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We report the first calculations of cross sections for the radiative capture reactions ^H(a,7)^Li 
and ■^He(Q,7)^Be below 2 MeV which use wave functions derived from realistic nucleon-nucleon 
interactions by the variational Monte Carlo technique. After examining several small corrections 
to the dominant El operator, we find energy dependences for the low-energy S-factors which agree 
reasonably with experimental measurements. There is no contradiction with the previous theoretical 
understanding of these processes, but the zero-energy derivative of the ^il{a, 7)^Li S-factor is smaller 
. than in most models. While this method can in principle predict cross section normalizations, the 

' normalizations of our results are mostly too low. 

(N 

■ I. INTRODUCTION 

[ Electromagnetic captures of alpha particles on and ^He are important processes in astrophysics. Together, they 
• are responsible for all ^Li production in the standard big bang nucleosynthesis (BBN) calculation. Because their 
^ I cross sections are also difficult to measure in the laboratory at the relevant energies (20-500 keV), they are the major 
I sources of uncertainty in the calculated primordial ""Li abundance '^He(Q;, 7) ''Be is also important for predicting 

psj ' the production rate of ^B and ''Be neutrinos in the sun. Accurate knowledge of its rate at solar energies (~ 20 keV) 
PsJ . is therefore important for studies of the solar neutrino problems ^ . 

' Recent theoretical studies of these cross sections have used two closely-related approaches. Potential models 

treat alpha particles and tritons (or '^He nuclei - the isospin symmetry of these systems makes the two processes almost 
identical) as point particles and the final states as bound states of these point particles. Authors of such models 
generate wave functions from potentials that fit experimental scattering phase shifts and bound state properties 
(binding energies, electromagnetic moments, etc.), and then compute the cross section as a direct capture. The 
1^ I resonating-group method (RGM) |p|,p|-p7t is fully microscopic, in that it solves an explicitly seven-body problem 
I , with a nucleon-nucleon potential, the parameters of which are adjusted to reproduce bound-state and resonance 
^ properties for the particular problem being solved. The name is derived from the choice of basis states for solving the 
^ , Schrodinger equation, which consists of one or more partitions of the particles into clusters with internal harmonic 
Ph • oscillator structure. The potential models now have a well-founded justification in the resonating-group work, in the 
j>! [ form of the microscopic potential model [|l8|,^ . Although the RGM models have apparently been very successful in 
describing these reactions (the calculation of Kajino correctly predicted both energy dependence and normalization 
for the experiment of Brune et al. |^), it is not clear that this is the final word, because the agreement with the data 
seems to be spoiled when the model space is expanded fl^,!!^ . 

It was shown in the early work of Christy and Duck ||2l|] that low-energy radiative captures on light nuclei may 
be treated to good approximation as external direct captures, that is, as one-step processes in which most of the 
matrix element arises outside the nuclear interaction radius. It remains true in more detailed models that the largest 
contributions to the matrix elements arise in regions well outside the range of the internuclear forces. In principle, 
then, the cross section energy dependence is given by convolutions of (positive- and negative-energy) Coulomb wave 
functions with the current operator, while the normalization of the cross section is determined by the asymptotic 
normalization of the bound state in the appropriate clusterization channel. 

However, fine details are obviously missing from the external direct capture model, motivating application of more 
detailed models that can predict both energy dependence (which is relatively easy to measure) and normalization 
(which is not easy to measure, and which must be inserted by hand in an external direct capture calculation after 
determination by some other means). Small effects involving the short-range (< 5 fm) behavior of the nuclear wave 
functions, both in the alpha-trinuclcon channel and in other channels, can affect the cross sections by several percent. 
Such effects are probably behind the differences between models pointed out in Ref. M^. In principle, there are also 
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small corrections to the leading-order current operators. Because of the astrophysical importance of determining these 
cross sections, and especially the need for low-energy extrapolation of '^He(a,7)^Be for solar physics, it is important 
to apply new approaches to this problem as they become available, and compare the results with past efforts. 

In this context, recent developments in the physics of light nuclei are particularly interesting. There now exist 
"realistic" nucleon-nucleon potentials which describe the np and pp scattering data, as well as the deuteron, with 
high precision {e.g. Ref. [|2^). Further interactions not describable by two-body potentials are described by three- 
nuclcon potentials, which have been adjusted to reproduce residual effects in the energy spectrum of light nuclei and 
properties of nuclear matter [p3| . Wave functions have been developed for these potentials in systems with up to 
eight nucleons . This provides an opportunity to approach the problem of radiative captures on light nuclei using 
realistic potentials and the computational techniques that have been developed to utilize them. 

Conversely, astrophysical interest in these processes has resulted in relatively precise measurements, which make 
cross section calculations useful tests of the wave functions in addition to the reproduction of the static moments, elec- 
tron scattering properties, and energy spectra to which they have been compared in the past. We have already reported 
the application of these wave functions to a radiative capture calculation in a paper on the process d{a, 7)^Li | p5| . 

The remainder of this paper describes cross section calculations for the reactions ^H(a,7)^Li and ^He(a,7yBe, 
using bound-state wave functions derived from realistic potentials by the variational Monte Carlo method. It is 
organized as follows: In Section ||, we describe the wave functions used to compute the cross sections. In Section 



[II, we describe the electromagnetic current operators and the methods used to compute their matrix elements. In 
Section [V, we describe the results for cross sections and branching ratios. In Section |V| we examine the implications 
of our results. 



II. WAVE FUNCTIONS 



A. Bound states 



We used ground states of ^H, '^He, ^He, ^Li, and ^Be which were found by the variational Monte Carlo (VMC) 
technique for the Argonne wis two-nucleon potential (hereafter AV18) and the Urbana IX three-nucleon potential 
(UIX) . The radiative captures can go to either the ground state or the (bound) first excited state in both ^Li and 
^Be, so the first excited states of these nuclei were also needed. These wave functions were generated by the same 
VMC method as the ground states. The bound-state wave functions have been reported in Refs. (triton and ^He) 
and (modified here as in Ref. to obtain ^Li and ^Be bound states with desired asymptotic properties). 

The VMC method proceeds by constructing wave functions as products of pair and triplet correlations between 
nucleons, and adjusting free parameters in these correlations to minimize energy expectation values which are com- 
puted by a Monte Carlo integration. The bound state wave functions are built from central and operator correlations 
between nucleons, acting on a Jastrow wave function, 

\^J) , (1) 

where Uij and Uf^^^ are two- and three-body correlation operators that include spin and isospin dependence and S is 
a symmetrization operator, needed because the Uij do not commute. The sums and products throughout this Section 
are over all nucleons. For '^He, the Jastrow part takes a relatively simple form: 

I*./) = n /y'' n /(^y) I*- (0000)1234) , (2) 

i<j<k<4 i<j<4 

where f{rij) and fijk are pair and triplet functions of relative position only, and $0(0000) is a determinant in the 
spin-isospin space of the four particles. Jastrow wave functions for ■^Hc and for the triton are constructed analogously, 
but the parameters of the /(r^) for these nuclei have been chosen to minimize energy expectation values of the 
three-body nuclei rather than of the alpha particle. The triton and '^He are identical in our calculation, except for 
their isospin vectors. In cases where the distinction is unimportant, we refer to both nuclei as "the trinucleon," and 
denote them both by r in subscripts that label clusters. 

For larger nuclei, spatial dependences must be introduced to place some particles in the p-shell. The A = 1 Jastrow 
wave function is constructed from scalar correlations multiplying a shell model wave function, 
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LS[n] 



where A is an antisymmetrization operator over all partitions of the seven particles into groups of four and three. 
For the central pair and triplet correlations, fxy{i'ij) and f^^f^\ the xyz denote whether the particles are in the s- 
or p-shell. The shell model wave function |$7(LS'[n] JAfTTJj)) has orbital angular momentum L, spin S, and spatial 
symmetry [n] coupled to total angular momentum J, projection A/, isospin T, and charge state T3, and is explicitly 
written as 



|$7(L5JAf[n]rr3)i234:567) 
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The Ylm{^) are spherical harmonics, xi^i'i^s) are spinors, and vlt^t^) are spinors in isospin, while brackets with 
subscripts denote angular momentum and isospin coupling. 

Particles 1-4 are placed in the s-shell core with only spin-isospin degrees of freedom, while particles 5-7 are placed 
in p-wave orbitals 0p (Rak) that are functions of the distance between the center of mass of the core and particle k. 
Different amplitudes PLS[n] in Eq. are mixed to obtain an optimal wave function; for the J'^, T = 3/2~, 1/2 ground 
state of ^Be, the p-shell can have /^iipj, /3ii[2i]i /3i|[2i]i /52f [2i]i and /32i[2i] terms. By far the largest contribution 
from these terms, as expected and as derived by diagonalization of the variational wave functions, is fi^i^^y This is 
true of all the mass-seven, T = 1/2 bound states, which are the final states of the radiative captures in question. 

The two-body correlation operator Uij is defined as: 



E 

P=2,6 



n f^3ki^^k.YJk) 



Up{r,j)Ol 



(5) 



Tj, and the /f^.^ 



is an operator-independent 



where the Of^- ~ Ti ■ Tj, (Ji ■ a-j, cTi ■ (TjTi ■ tj, Sij, and Sij'i i ■ • j, 

three-body correlation. The six radial functions fssif) and Mb=2.6(^) are obtained from two-body Euler-Lagrange 
equations with variational parameters as discussed in detail in Ref. |2q| . They are taken to be the same in the p- 
shell nuclei as in "^He, except that the Up=2,6{'i^) are forced to go to zero at large distance by multiplying in a cutoff 



factor, 



1 + exp[-i?„/au] / 1 + exp[(r - i?„)/at, 



with Ru and variational parameters. The fsp correlation is 

constructed to be similar to fss for small separations, but goes smoothly to a constant of order unity at large distances 
(r > 5 fm): 



fsp{r) 



1 + exp[(r - Rsp)/asp] 



fssir) + Cspil - exp[-(774p)^]) 



(6) 



where asp, b^p, etc. are additional variational parameters. The fpp{r) correlation in the mass-seven nuclei is the same 
as the f(rij) correlation in the trinucleon, so that when the three p-shell nucleons are all far from the s-shell core, 
they look very much like a trinucleon. 

These choices for /^s, fsp, fpp, and Up=2.6 guarantee that when the three p-shell particles are all far from the s-shell 
core, the overall wave function factorizes as: 



^"7 



(7) 
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where tpa is the variational "^He wave function, ipr is the variational trinucleon wave function, and Tar denotes the 
separation between the centers of mass of the a and trinucleon clusters. Provided that [fspirar)]^'^ goes to a constant 
quickly enough and smoothly enough, the long-range correlation between clusters is proportional to [4'p^^^\raT)]'^ ■ 

The single-particle functions (j)p^^^Rak) describe correlations between the s-shell core and the p-shell nuclcons, and 
have been taken in previous work [ p7| to be solutions of a radial Schrodinger equation for a Woods-Saxon potential 
and unit angular momentum, with energy and Woods-Saxon parameters determined variationally. It is important 
for low-energy radiative captures that these functions reproduce faithfully the large-separation behavior of the wave 
function, because the matrix elements receive large contributions at cluster separations greater than 10 fm. In fact, at 
20 keV, more than 10% of the cross section for •^He(Q;, 7)''Be comes from cluster separations beyond 20 fm. At these 
distances, well outside the nuclear interaction distance, the clusterization with the lowest cluster-separation energy 
should be the most important. We have therefore modified the bound-state wave functions for the capture calculation 
to enforce cluster-like behavior, matching laboratory cluster separation energies, when the three p-shell nucleons are 
all far from the s-shell core. 

In general, for light p-shell nuclei with an asymptotic two-cluster structure, such as ad in ^Li or at in ''Li, we want 
the large separation behavior to be 



(8) 



where Wkm{'2'yr) is the Whittaker function for bound-state wave functions in a Coulomb potential (see below) and n 
is the number of p-shell nucleons. We achieve this by solving the equation 




ii£ + l) 



+ V{r)+A{r) 



r0^^["l(r) = O, 



(9) 



with £ = 1, /i4i the reduced mass of one nucleon against four, and V{r) a parametrized Woods-Saxon potential plus 
Coulomb term: 



V{r) = - 



Vo 



2{Z-2) 



- cxp[(r - i?o)/ao] 



F{r) . 



(10) 



Here Vq, i?o, and uq are variational parameters, (Z — 2)/n is the average charge of a p-shell nucleon, and F{r) is a 
form factor obtained by folding a and proton charge distributions together. The A(r) is a Lagrange multiplier that 
enforces the asymptotic behavior at large r, but is cut off at small r by means of a variational parameter cq: 



A{r) = X{r) 1 -exp( - (r/co) 



(11) 



The A(r) is given by 



A(r) 



2fi 



41 



1 (fuL 2 



UL 



2{Z-2) e2 



(12) 



where is directly related to the Whittaker function (solution to the radial Schrodinger equation for negative-energy 
states in a purely Coulomb potential): 



WL/r= (T^fc™(27r)/r)i/'^ 



(13) 



Here 7^ = 2/i4„_B4„/fi,^, with /i4„ and B^n the appropriate two-cluster reduced mass and binding energy, k = —2[Z — 
2)e^/i4„//i^7, and m = L + i. 

For the ''Li ground state, the largest contribution has ^43 = 2.47 MeV (binding energy of ''Li relative to a and 
t clusters) and L = 1 corresponding to the asymptotic P-wave of the ''Li ground state, or amplitude /^iipj in Eq. 
(^. None of the other possible amplitudes /3LS[n] correspond to asymptotic at clusterizations. However, there is no 
reason for them not to be present in compact configurations of the nucleons. Including such components in the wave 
functions improves the binding energies of the mass-seven bound states by about 0.2 MeV. The asymptotic forms of 
^LS[n] lower-symmetry channels are set to match the threshold for ''Li ^Li -I- n. Analogous descriptions 

hold for the other bound states (the ''Li excited state and the two ''Be bound states) , with the appropriate thresholds 
substituted for ^43. 
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The ''Be ground and first excited state Jastrow functions have been treated in previous development of the variational 
Monte Carlo wave functions |^ as the isospin rotations (T3 = +1/2 instead of —1/2) of the corresponding ''Li shell- 
model-like wave functions. In this work, the ''Li and 'Be bound state Jastrow functions also differ by the choices 
of Bin for the asymptotic cluster behavior of the wave functions, which match the cluster breakup thresholds as 
described above in each case. This choice of B^n has the result that in configurations where the p-shcU nuclcons 
are far from the s-shell core, the energy is the sum of the Coulomb potential, the kinetic energy contributed by the 
0p "^'"^ , and the cluster binding energies (well-reproduced because the core resembles an alpha particle and the p-shell 
has been constructed to resemble a trinucleon). Because the t^p"^'"' matches the laboratory binding energy for the 
known Coulomb potential, the local energies at large particle separations match the known binding energies in the 
mass-seven nuclei. This agreement has been confirmed numerically. 

In Refs. [^,^ the /^y^f three-body correlation of Eq.(H) was a valuable and computationally inexpensive improve- 
ment to the trial function, but no f^j^ or f^^^ correlations could be found that were of any benefit. However, for the 
types of wave functions used here, it is found that the correlation 

= 1 + qi[fssinj)/fppin,) ~ l]exp[-g2(r„. r„,)] , (14) 

with i, j labels of p-shcU nucleons and gi,2 as variational parameters, is very useful ||2^. It effectively alters the central 
pair correlations between pairs of p-shell nucleons from their trinucleon-like forms to be more like the pair correlations 
within the s-shell when the two particles are close to the core. This correlation improves the binding energy by ~ 0.25 
MeV in ''Li. 

The authors of Refs. [ p7|j29| reported energies both for the trial function ^'t of Eq.(|l]), and for more sophisticated 
variational wave functions, which add two-body spin-orbit and three-body spin- and isospin-dependent correlation 
operators. The ^'y gives improved binding compared to ^'t in both the mass-four and the mass-seven wave functions 
considered, but is significantly more expensive to construct because of the numerical derivatives required for the 
spin-orbit correlations. In the case of an energy calculation, the derivatives are also needed for the evaluation of 
X-dependent terms in AV18, so the cost is only a factor of two in computation. However, for the evaluation of other 
expectation values the relative cost increase is ~ Thus in the present work we choose to use for non-energy 
evaluations; this proved quite adequate in studies of ^Li form factors |3C| ] and of the six-body radiative capture 
d(a,7)^Li H. 

The variational Monte Carlo (VMC) energies and point proton RMS radii obtained with are shown in Table || 
along with the results of essentially exact Green's function Monte Carlo (GFMC) calculations |^,^,|l| and the 
experimental values. We note that the underbinding of the A = 1 nuclei in the GFMC calculation arises from the 
AV18/UIX model and not the many-body method; it can be improved by the introduction of more sophisticated 
three-nucleon potentials . 

Although the present variational trial functions with the imposed Coulomb asymptotic correlations produced a 
variational improvement in the case of ^Li, they give approximately the same energies as the older shcU-model-like 
correlations of Refs. (2^j2^ for 'Li and 'Be. Unfortunately, because the variational energies of the A — 1 bound 
states are not below those of separated alpha and trinucleon clusters, it is possible to lower the energy significantly by 
making the wave functions more diffuse. Therefore, the variational parameters were constrained to give RMS charge 
radii which agree reasonably with experiment, as seen in Table |. Given (/ip'^^"', the easiest way to constrain the RMS 
radius in the variational procedure is to choose the parameters of the /sp correlations between s- and p-shell nucleons 
to adjust the probability of finding the p-shell nucleons far from the s-shell core. There was considerable freedom in 
the specific form of these correlations as long as asymptotic properties of the wave functions were not being tested, 
because of the insensitivity of energy expectation values to the tails of the wave functions. However, the form of the 
correlation in the ar channel depends on /j^ , as seen in Eq. (Q), so the large-cluster-separation parts of the wave 
function are very sensitive to the choice of /sp. (The twelfth power arises because there are four particles in the 
s-shell and three in the p-shell, and thus 4 x 3 = 12 sp pairs.) Prior to this work, VMC wave functions had cluster 
distributions that dropped by an extra factor of two beyond 5 fm, relative to the drop expected on the basis of the 
clusterization arguments given above. The extra drop had two closely-related effects: the asymptotic normalization 
coefficients (see below) for two-cluster breakup were too small, and so were the cross sections which we initially 
computed from them - by a factor of about two. The present wave functions perform more poorly than previous 
VMC wave functions on the ordering of the mass-seven bound states (a perennial difficulty for both VMC and GFMC 
because of the close spacing of the states), but they give larger asymptotic normalizations at a reasonable cost in 
binding energy. The relationship between nuclear size (as measured by quadrupole moments) and cross sections for 
direct radiative captures has been noticed before, and applied usefully both to these |Q and other reactions, most 
notably ''Be(p, 7)^B [ p3[. We note that our S'(O) and quadrupole moments for both mass-seven systems fit the general 
trends shown in Ref . Ill 711 . 
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TABLE I. Calculated VMC, GFMC, and experimentally measured energies, point proton RMS radii, and quadrupole mo- 
ments of ^H, '^He ''He and bound states of ^Li and '^Be. Numbers in parentheses are Monte Carlo statistical errors. 



Nucleus 


Observable 


VMC *T 


GFMC 


Experiment 


xl 


E 


— 0.10(^1 ) 


— 0.4/ 1^1 J 


— o.4o 






1.60(1) 


1.59(1) 


1.60 


xie 


Hi 






7 79 










1 77 


*He 


E 


-26.89(3) 


-28.34(4) 


-28.30 






1 /I O / 1 \ 

i.4o(i) 


1 /I c / 1 \ 
i.4o(i) 


1 /I o 


^Li 


E 


-31.26(8) 


-37.78(14) 


-39.24 






2.30(1) 


2.33(1) 


2.27 






-3.7(2) 


-4.5(2) 


-4.U6 




E 


-29.55(8) 


-36.23(14) 


-37.60 




/ 2\l/2 


2.41(1) 


2.52(1) 






O 


-5.9(3) 


-7 5('3") 




''Li* 


E 


-31.37(8) 


-37.53(15) 


-38.77 






2.35(2) 


2.35(2) 




^Be* 


E 


-29.70(8) 


-36.01(15) 


-37.17 






2.46(2) 


2.54(2) 





TABLE IL Asymptotic normalization coefficients (in fm"'^^^) for the overlap between mass-seven bound states and ar 
clusterization, computed from the VMC wave functions. Best available estimates are presented for comparison on the lower 
line. 





^Li 


^Li* 


^Be 


^Be* 


VMC 


3.4 ±0.1 


2.65 ±0.10 


3.55 ±0.15 


2.9 ± 0.1 


literature 


3.55 ±0.27 


3.14 


4.79 


4.03 



The asymptotic two-cluster behavior of the seven-body wave functions can be studied by computing the two- 
cluster ar distribution functions, {Aipa'4'™\Vat \ V'lT) ^^i, and its analogs for the other seven-body bound 
states we consider. These functions are described in Ref. [Q. They can be expressed in terms of Clebsch-Gordan 
factors, spherical harmonics, and the radial functions Rijar) plotted in Figs. |^ and ^ At large Tq,-, rarR{i'aT) 
should be proportional to a Whittaker function, as described above. The proportionality constant is the asymptotic 
normalization constant Ci. We have extracted Ci from the overlap functions by a least-squares procedure, matching 
them to the appropriate Whittaker functions. We find that for all the mass-seven wave functions considered here, 
Ci becomes asymptotic at r^T ~ 7-9 fm. In Table |l| we present asymptotic normalization coefficients for the A = 7 
bound states, based on fitting overlaps in the region 7-10 fm as shown in Figs. |l| and |[ 

Our values of Ci for the mass-seven nuclei tend to be somewhat smaller than the values found in the literature . 
Only in the case of the ^Li ground state are experimental determinations of the Ci of reasonable quality. Brune et al. 
p5| find a "world average" of 3.55 ± 0.27 fm~^/^, relying mainly on theoretical models djlj], and giving less weight 
to the partially experimentally-based evaluations of Igamov |36|. For the other states, we rely on the Igamov et al. 
p6| extraction of ANCs from Kajino's RGM calculations with the MHN potential These calculations match the 
radiative capture data for ■'H(Q:,7)^Li very well in both energy dependence and normalization. However, Igamov et 
al. report "nuclear vertex constants," which differ from ANCs by prefactors whose definitions are ambiguous in the 
literature; the numbers presented in the second row of Table || should be used with caution. The present results for 
the ANCs are also subject to correlated uncertainties characteristic of the Monte Carlo integration method used to 
compute the overlaps, and the uncertainties ar e the refore difficult to estimate reliably (as discussed in more detail 
with regard to matrix element densities in Sec. Ill A below). 

As noted in the Introduction, low-energy direct captures can be treated to good approximation by considering 
only cluster separations beyond a few fm, and using only the longer-range parts of the bound- and initial-state wave 
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FIG. 1. Monte Carlo samples of the radial two-cluster at distribution functions in ^Li ground state [J^ = 3/2~, left panel) 
and first excited state (J"^ = l/2~, right panel), with error estimates. The solid curve is the expected Whittaker-function 
asymptotic form of the overlap, normalized as in Table |l| to match the ar distributions at 7-10 fm. 
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FIG. 2. Same as Fig. |l| but for radial two-cluster a^He distribution functions in ^Be bound states. 

functions to compute matrix elements. This approach requires the provision of "spectroscopic factors," or more 
precisely, the ANCs discussed above. Because we have imposed the condition that the large-cluster-separation part 
of the ground state match its expected form, it is true that for our cross section calculations, the VMC method 
provides Ci. However, it also provides the inner few fm of the wave function, which requires some model of what is 
going on inside the nuclear interaction radius, and may be important for understanding differences in the logarithmic 
derivatives of the S'-factor found in various theoretical studies. 





B. Scattering states 

The initial-state wave functions are taken to be elastic-scattering states of the form 

\4>^r;LSJM)=AUi^ir^r)YLMA^c.r)]jG,,\^ba.4'T')\ , (15) 

I ) LSJM 

where curly braces indicate angular momentum coupling, A antisymmetrizes between clusters, V'a is the *He ground 
state, and iplf}^ is the trinucleon ground state in spin orientation ms- 

The Gij are identity operators if the nucleons i and j are in the same cluster. Otherwise, they are a set of both 
central and non-central pair correlation operators which introduce distortions in each cluster, under the influence of 
individual nucleons from the other cluster. They are derived from solutions for nucleon-nucleon correlations in nuclear 
matter , and become the identity operator at pair separations beyond about 2 fm. (These correlations have been 
included in the definition of the overlap functions shown in Figs. ^ and ^.) 

The correlations are derived phenomenologically. The variational seven-body bound-state wave functions do 
not give the correct energies with respect to cluster breakup, so we do not expect to be able to use the variational 
technique to solve for these correlations. Instead, we generate the (f)'l^^ as solutions to Schrodinger equations, from 
cluster-cluster potentials that describe phase shifts of a-^H and a-'^He scattering as scattering of point particles. 
Because of the small amount of available laboratory data and the large amount of work that has already been put 
into generating potentials that reproduce them, we take cluster-cluster potentials from the literature p|-[^,pj|]. We 
now point out the main features of these potentials. 

We treat all of these models as (and many have been explicitly constructed as) descriptions of both the at and 
a ^He systems, with appropriate Coulomb potentials and laboratory masses. Each of the potentials we use to generate 
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the (pi^^ has a deep, attractive central term, and a spin-orbit term. The spin-orbit terms are constrained mainly by 
the spacing between the P3/2 and P1/2 bound states, and between the resonances at 2.16 and 4.21 MeV in F7/2 and 
F5/2 at scattering, respectively (as well as their analogs in a'^He scattering). The other interesting feature in the 
scattering of the odd partial waves is the apparent hard-sphere behavior of the phase shifts in P-wave scattering. 
This comes about because the wave functions between the clusters must respect the Pauli principle by allowing 
antisymmetrization of nucleon wave functions between clusters. In the P waves, this takes the form of wave functions 
that have a single node whose location is almost independent of scattering energy, and it therefore gives rise to phase 
shifts which look like scattering from a hard sphere |Q. (See Figs. and g, where the steep dips in absolute values of 
the cluster distributions at 2 fm correspond to nodes). In the cluster-cluster potentials, this requires that the central 
term be large and negative, so that the ground state wave function has one node. Of course, a more tightly-bound 
"forbidden" state with zero nodes exists for such a potential, but it does not allow antisymmetrization of the nucleon 
wave functions. We note that the requirement of a particular nodal structure only models approximately the effects of 
the Pauli principle on the inter-cluster correlations. We did not use potentials that enforced the hard-sphere behavior 
with repulsive short-range terms |39| . 

Because of the requirements of the Pauli principle, one also expects different potentials to describe the odd- and 
even-parity scattering. The even-parity phase shifts are unfortunately lacking in details which models must match 
(see Fig. ||), beyond the apparent hard-sphere scattering in the 5*- wave data, corresponding to the Pauli-required 
minimum of two nodes in the wave function [Q. The measured D-wave phase shifts are even worse, being consistent 
with zero (or 180°, with the knowledge that there must be at least one node in the wave function) throughout the 
region below 8 MeV for both a-t and a-^He systems pO[-^. This lack of features is particularly unfortunate in that 
by far the most important reaction mechanisms for radiative capture in this system are El captures from S- and 
ZJ-wave scattering states. The most stringent tests of any of the D-wave potentials have been comparison of their 
phase shifts with the results of more elaborate (RGM) theoretical models. 

Finally, we note the difficulty of reproducing the published models, which is due to omitted descriptions of details 
of the potentials, particularly handling of the Coulomb potential at short range. Rather than guess how to fix up each 
potential, we restrict attention to potentials that allowed good descriptions of the low-energy scattering on the first try 
(with no short-range cutoff in the l/r Coulomb term unless explicitly described by the potential's authors). Relatively 
small differences in the cluster-cluster potentials are directly connected to the size of the scattering wave function at 
energies less than 1 MeV and cluster separations less than 20 fm, resulting directly in differences in the normalization 
of our computed capture cross sections from one potential to the next. However, it is possible to eliminate the worst 
potentials on the basis of the low-energy phase shifts; potentials which underpredict the phase shifts also produce 
radiative capture cross sections that are too low by as much as a factor of two, relative to those generated from other 
cluster-cluster potentials. Potentials which were created for use in orthogonal cluster models (as opposed to simple 
cluster models) were eliminated from application to our problem on this criterion. 



III. OPERATORS 

The cross sections were computed by a multipole expansion of the electromagnetic current operator H4 



E 



LSJiJ;. 



2Jf + 1 Vrei l + q/mr 



M. 



LSJ^Jf 



(9) 



(16) 



where a is the fine structure constant (a ~ e^/hc), Wj-ci is the ar relative velocity, my is the mass of the final state, 
and Ef^'^^'^^{q) and M^^'^^'^' {q) are the reduced matrix elements (RMEs) of the electric and magnetic multipole 
operators with multipolarity i connecting the scattering states in channel LSJi to bound states of ''Li or ^Be with 
angular momentum Jf. The center-of-mass energy of the emitted photon is given by 



my 



1 + W IH (m^ + TUa - in-! + Ec.m.) 

V TTij 

llr - mj + i?c.m. 



(17) 



where m^, m^, and my are the rest masses of the trinucleon, ^He, and the appropriate ^Li or ^Be state, respectively. 
The astrophysical S'-factor is then related to the cross section via 



S{Ec.m.) = i?c.m. Cr(i?c.m.) exp(2ZiZ2 Tra/ftUrol) , 



(18) 
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FIG. 3. Phase shifts produced by the potentials used to generate the cjiirr- Data are taken from Refs 
panel, x denotes J'^ = |^ phase shifts, o denotes J'^ = I""" phase shifts. Potentials are taken from Refs. 
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where Zi and Z2 are the charges of the two initial-state nuclei. 

The dominant reaction mechanism, for creation of both the excited state and the ground state, is an El (electric 
dipole) transition from an S-wave scattering state to a bound state. Capture into the excited state is followed 
immediately by electromagnetic decay to the ground state, so the total cross section of interest for astrophysics is the 
sum of the cross sections for captures into the two states. At energies above about 500 keV, E\ capture from D waves 
becomes important. We computed transitions originating from scattering states with orbital angular momentum 
L = 0,1,2 and 3, via Ml, M2, El, E2, and E3 transitions, and found that up to the 0.1% level, only El captures 
originating from S- and D-wave scattering states matter at energies below 1.5 MeV. 

With the exception of the El term, all RMEs were computed in the standard long- wavelength approximation (LWA), 
keeping only the lowest-order term in photon wavenumber of the modified spherical Bessel functions appearing in the 
RME integrals. The LWA is valid to reasonable accuracy because at the low energies under consideration {E < 1 
MeV), the ratio of system size to photon wavelength is less than 10 fm/200 fm = 0.05. 

In a previous study p^], we developed code to examine the isospin-forbidden El transition in d(Q;,7)^Li. We have 
applied this code to compute corrections to the LWA for the El transitions under consideration here. An examination 
of these corrections for the mass-seven system is in principle of interest for the problem of extrapolating cross sections 
to low energies. However, we find that all but one of them provide contributions of less than 0.05% of the total cross 
section. (See Ref. |^ for a list and detailed discussion of the corrections we applied, which extend to third order in q.) 
We do not actually compute the largest correction, which is the "center-of-energy" correction. This correction arises 
because potentials and kinetic energies should be included in the definition of the center-of-momentum frame, but have 
not been The center-of-energy correction becomes important when the leading order LWA operator vanishes, 

as in c?(q!, 7)^Li, but it should amount to only +2.4% of the El cross section in ■^H(a,7)^Li and -1-3.1% of the El 
cross section in ■^He(Q;, 7) ''Be. We did not compute center-of-energy corrections (or include an estimate of them in the 
results presented below) because of the extra computation necessary to find energies during the capture calculation. 
Their omission is not serious because 1) their effect is to change the normalization, not the energy dependence of the 
El cross section, 2) our model calculation is only accurate to about 5-10% at best, and 3) the above estimates of the 
size of the effect should be quite accurate, being based only on the differences between using nuclear masses and using 
integer multiples of the mean nucleon mass in the factor [(Z2TO1 — Zim2) / [mi + m2)Y in the LWA cross section. 



9 



A. Matrix element integration 



Actual computation of the matrix elements was performed with a modified version of the code described in Ref. P5[ , 
which is itself a modified version of a code developed to compute energies and other properties of light nuclei |27| ] 
for variational calculations. The method used to integrate over nucleon configurations is the Metropolis Monte Carlo 
algorithm, with a weight function proportional to the bound state wave function involved in the computed transition. 
As discussed below, this weight function was chosen to reflect in general detail the form of the matrix element 
integrands, and to obtain significant numbers of Monte Carlo samples over a broad range of cluster separations. The 
final calculation consisted of 10^ samples for each transition. 

We have applied the approach of splitting the calculation into energy-dependent and energy- independent parts, as 
in Ref. , so that the reduced matrix elements of Eq. (^6|) are written 



' <J,M,,eX\Jfmf> ^ ''> 

X dxx^ (/.fv^(x)(V'7^'"^ mx{q)A I S{X - r^r)Y^'' (far)) J] Gyl^aV^r" ) > (20) 

I ) LSJiMi 

and computed using photon polarization A = -f 1 for the multipole operators Te\{q). The delta function is applied by 
accumulating the Monte Carlo integral in radial bins of thickness 0.25 fm. The final integration over x is performed 
by inserting the appropriate dependences on photon energy in each term (since the LWA expansion is in powers of 
energy, this dependence may be taken out of the integral), and computing (j)"^!^, at each energy. This allows the time- 
consuming Monte Carlo integration to be performed only once for each partial wave and operator, so that computation 
of RMEs for many energies is relatively inexpensive in computer time. After initially setting up the code and checking 
that selection rules were satisfied, we did not explicitly compute RMEs for parity-forbidden operators. 

At cluster separations beyond about 10 fm, the RME densities were subject to considerable noise in the Monte Carlo 
sampling. This is because while Monte Carlo weighting schemes based on the ground state give good sampling along 
directions other than the cluster separation in the 3A-dimensional configuration space, they provide small numbers 
of samples at large Tqt- (due to exponential decay of the wave function at large distances) . VMC work usually uses 
weighting based on the inner product in spin-isospin space of a simplified bound state with itself. Such weighting is 
good for minimizing Monte Carlo variance of integrands that resemble the bound state probability density closely, such 
as energies, but it does not provide enough samples in the asymptotic region of the wave function to do low-energy 
direct capture calculations. Our previous paper ||2^ ] used the square root of the usual weight function, extending 
the Monte Carlo sampling out to large cluster separations, but at the expense of greater sampling noise for a fixed 
number of samples at fixed cluster separation. In the mass-six problem, it was straightforward to run the Monte 
Carlo integration of the RMEs until the densities had "converged" to the expected asymptotic forms at large cluster 
separation. This required about 2 x 10^ total configurations for a given RME. At ^ = 7, the spin-isospin space is 
larger, so the code is slower by a factor of about 10, and it was only practical to obtain 10^ samples for each RME 
with available computing resources. However, since the configuration space gains three dimensions with the addition 
of a particle, 10^ samples do not provide as thorough of a sample of the configuration space at mass 7 as at mass 6. 
The result is at best a few-percent measure of the asymptotic normalization at 7-10 fm. Many samples are obtained 
at larger cluster separation with the new weighting scheme, but the samples beyond 15 fm appear to have correlated 
noise. This is presumably because the samples in question correspond to only a few excursions of the sampling Markov 
chain into the region of large cluster separation, and are therefore not very independent from each other. They tend 
to be either mostly high or mostly low relative to the asymptotic forms explicitly built into the wave functions. See 
Fig. U for illustration of these difficulties. A much larger sample (by a factor of 10) would be expected to exhibit 
much less of this sort of correlated noise, but is prohibited by the large amount of computer time that would require. 

The matrix elements were therefore integrated out to 7 fm using the Monte Carlo results for the integrand of Eq. 
(pj|); integration beyond 7 fm was carried out using the known asymptotic forms of the matrix element densities, 
normalized to the Monte Carlo output by least-squares fitting at 7-10 fm. This range was arrived at by comparing 
results for the ANC of the ^Li ground state arising from two different weighting schemes and varying numbers of 
samples. Reasonably consistent agreement was found by fitting at 7-10 fm in all cases. The accuracy in the cross 
section then depends on how accurately the ANCs can be determined in this region. 



IV. CROSS SECTIONS 
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FIG. 4. Ratios of the computed radial overlap functions of Figs. |l| and I to their imposed asymptotic forms, based on 
lO" particle configurations. These ratios should be equal to the asymptotic normalization constant for the appropriate wave 
function beyond about 7 fm. Superimposed on the ratios for the ^Li excited state is the corresponding distribution of Monte 
Carlo samples, essentially identical for all four overlap functions. 
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FIG. 5. Low-energy S-factor for ■^H(Q:,7)^Li (thick solid line), and its breakdown into contributions from various partial 
strengths, computed using potential A of Kim et al. [^. Transitions to the ground state are shown as solid lines, and transitions 
to the excited state as dashed lines. Labels indicate the initial state, and additional symbols indicate multipole operator: no 
symbol. El; o, E2; x, Ml; □, M2; and A, order spin correction to El operator (LWA2 of Ref. [|||). 



A. ='H(a,7)^Li 

The computed S'-factor for '^H(a,7)''Li is shown broken down into contributions from various terms of Eq. ( p^ in 
Fig. and in comparison with laboratory data in Fig. |6[ The dominant processes are obviously El captures, with 
large contributions from captures into both the ground and excited states. Contributions from other partial waves 
and multipole operators are not present above the 1% level. Contributions from higher-order LWA corrections to the 
El operator are less than 0.02%. Our calculations are therefore limited in accuracy only by 1) the accuracy of the 
(bound- and scattering-state) wave functions, and 2) Monte Carlo statistics. 

The present calculation is the result of sampling 10^ points in the seven-particle configuration space. Formal error 
estimation on the resulting cross sections is difficult because all 35 partitions of the seven nucleons into alpha and 
triton clusters were computed for each configuration. This saves computation time, and enforces antisymmetry of 
the initial state exactly; it also introduces correlations between values of the operator densities (integrands of Eq. 
(|l9|)) at different cluster separations. The uncertainties in different transitions and partial waves are also correlated, 
because they are based on the same random walk of particle configurations. There is also an uncertainty from the 
wave function normalizations, because they are also derived from Monte Carlo integrations. This last uncertainty 
amounts to about 3% in the cross section. 

We take the best indication of the Monte Carlo uncertainty to be the formal uncertainties on the asymptotic nor- 
malizations of the Monte Carlo matrix element densities, which we used to compute matrix element contributions 
beyond cluster separations of 7 fm. The asymptotic forms of the matrix element densities are the Whittaker-function 
asymptotic forms discussed above, multiplied by the radial dependences of the electromagnetic multipole operators. 
The normalizations were fitted to the matrix element densities at cluster separations of 7-10 fm. Using fitted asymp- 
totic forms amounts to treating a large part of the matrix element as an external direct capture, and it fixes two 
problems. First, it removes the need for large numbers of Monte Carlo samples in the remote tails of the wave function. 
Second, the asymptotic normalizations are found by a weighted least-squares procedure, and formal error estimates 
on these normalizations are possible. In practice, the correlations between matrix element densities produce reduced 
significantly less than unity. A common approach when confronted with such a problem in experimental data is to 
assume that uncertainties have been overestimated, and to reduce formal uncertainties accordingly. We have not done 
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FIG. 6. Same as Fig. ^, but showing only total S'-factors for several different potentials used to generate cluster correlation 
functions : dot-dashed and dotted, potentials A and B respectively of Kim et al. j^; short-dashed, from Dubovichenko et 
al. solid, from Buck et al. and long-dashed, from Buck and Merchant [^. Data are from Brune et al. [^, and share a 
common 6% normalization uncertainty, not shown. 



this, and we arrive at uncertainties of approximately 10% in the cross section, based on the formal error estimates 
(corresponding closely to the sizes of ANC errors in Table More detailed analysis is problematic, and it is not 
called for because the 10% estimate is already larger than any other contribution to the error budget. 
The results themselves are best characterized in three ways: 

1) Normalization. Our results are lower than the data of Brune et al. |Q by 0-20%, depending on the cluster 
potential. Although other data exist p7|-pT| , those of Brune et al. are much more precise, and permit the best test of 
our results. Those data share a common 6% normalization uncertainty, not shown in Fig^. The systematic discrepancy 
between some of our results and the data is probably small enough to ascribe to the combined uncertainties of the 
Monte Carlo integration and of the data. Taking the branching ratio to have its experimental value of i? = 0.453, 
the computed S'-factors for transitions to the ground state match the measured low-energy S'-factors, while those for 
the excited state do not. The normalization is affected significantly by the choice of potential used to generate the 
inter-cluster correlations for the scattering states. By applying five different potentials from the literature as 
described above, we obtain a variation of ±5% in the 5(0) (total range for the five potentials) about a mean of 0.90 
keV-b - a full range equal in size to the Monte Carlo uncertainty. Summarizing our results in one number, we obtain 
using the inter-cluster potential A of Kim et al. Q (the best fit to the S'-wave phase shifts), 5(0) — 0.095 keV b, 
similar to other estimates found in the literature. 

2) Branching ratio. The branching ratio i?, defined as the ratio of the cross section for capture into the excited 
state to that for the ground state, is shown in comparison with the Brune data in Fig. 0. A weighted least-squares 
normalization of the calculation to the laboratory data shows that our calculation of R is lower than the data by 15%, 
which can be combined with the results above to infer that the excited-state ANC is low by 8%, within the range of 
the Monte Carlo sampling error. The energy dependence of the branching ratio matches the data with a of 26.0 
for 15 degrees of freedom, about as well as the straight line fits of Brune et al. 

3) Energy dependence. The energy dependence of the ■^H(a,7)^Li S'-factor at low energy is almost completely 
independent of cluster-cluster potential, for the five potentials examined. After normalizing the computed cross 
section (with cluster correlations computed from potential A of Kim et al. [^) to match the Brune data, we obtain a 

of 38.7 for 16 degrees of freedom. A significantly better fit results if the highest 3 points in energy (where D-wave 
capture becomes important) are excluded. The residuals of —1.8% to +8.5% compare favorably with other theoretical 
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FIG. 7. Computed ratio of the cross section for capture into the excited state to that for capture into the ground state, for 
^H(a,7)^Li. Symbols are the same as in Fig. ^, and the data share a normalization uncertainty of 4% (not shown). 



calculations, but the energy dependence of the calculation is systematically shallower than that of the data at the 
highest-energy points. Our calculation of the S factor gives a logarithmic derivative at = of —0.972 MeV^^ - 
about equal to that found by Mohr et al. but half of that found in other theoretical studies |Q,|l^,0, and about 
2/3 of that suggested in a recent compilation of astrophysical reaction rates 



B. ^He(a,7)''Be 

The low-energy S'-factor computed from the Kim A cluster-cluster potential for •^He(a,7)^Be is shown in Fig. ||, 
along with the contributions of individual terms of Eq. (|l^). Total S*- factors from five cluster-cluster potentials 
are shown along with the laboratory data in Fig. ^ The computed branching ratios are shown along with the 
corresponding data in Fig. |l^. The discussions of small contributions to the S'-factor and of the precision of the 
results for ■^H(a,7)^Li above also apply here, again with a Monte Carlo error estimate of about 10% in S'-factor 
normalization. We again break down the results for ■^He(Q;, 7)''Be into normalization, branching ratio, and energy 
dependence: 

1) Normalization. Our S'-factors are more than 10% lower than the lowest data set (after applying the renormaliza- 
tion of the Krawinkel data set recommended by Hilgemeier et al. and nearly a factor of two lower than the 

highest data sets (note that this includes the data of Volk et al. , which are not shown because their results were 
published only as extrapolated S'(O) values). By applying five different phenomenological cluster-cluster potentials 
that are not in disagreement with the low-energy elastic-scattering data, we obtain a small variation in S(0) about 
a mean of S'(O) = 0.40 keV b. Using the potential which best matches the low-energy S'-wave scattering (potential 
A of Kim et al. ||^), we obtain S'(O) = 0.40 keV b. Although this is closer to matching the lower numbers found in 
capture photon experiments than the delayed activity experiments, it is not a close match in normalization to any of 
the experimental results. Our results are therefore not useful for addressing possible systematic problems in the data. 



2) Branching ratio. It is seen from the branching ratios in Fig. 10 that our calculation is in reasonable agreement 
with the laboratory data with regard to relative strengths of transitions to the ground and first excited states of ^Be. 
This suggests that the ANCs of the VMC wave functions for both states are too small, and by about the same factor 
in each case. Using the results for S-factor normalization above, we conclude that this factor is in the range 5-25%. 

3) Energy dependence. We renormalized our results to best fit each of the larger data sets separately, and computed 



chi-squared statistics in each case to determine goodness of fit. The results are shown in Table [II, and indicate general 
agreement. For the logarithmic derivative of S{E) at E = 0, we obtain —0.57 MeV~^, in reasonable agreement with 
other estimates in the literature ||7|, |l4| , pT|j52[ , where published models fall in the range —0.50 to —0.72. It is about 
equal to the value presently used in solar neutrino work |^]. 
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FIG. 9. Same as Fig. ^, but showing only total S'-factors for several different potentials used to generate cluster correlation 
functions 4>ai, with symbols as in Fig. ^ Data are from Refs. [p3|-p8[. Symbols for the data are the same as in Ref. with 
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FIG. 10. Computed ratio of the cross section for capture into the excited state to that for capture into the ground state, for 



^He(a, 7)^^66. Symbols are the same as in Fig. ^ 



TABLE III. Comparison of energy dependence of calculated 5-factors with laboratory data, after renormalizing the computed 
S- factor to minimize x^- v denotes number of degrees of freedom, and an estimate of systematic normalization uncertainty was 
subtracted from the error estimates of the points (as discussed in Ref. ||l]) before performing this analysis. 

Data set ~)i 1'^ ^ 

Krawinkcl @ 0.364 37 

Parker Q 0.917 37 

Hilgemeier jEl 0.698 8 

Nagatani Ej 0.450 6 



16 



V. IMPLICATIONS AND APPLICATIONS 



We have carried out the first calculation of the low-energy S'-factors for the processes '^H(q;, 7)^Li and ■^He(a, 7)^Be 
based on realistic nucleon-nucleon potentials. Seven-body wave functions for these potentials, constructed by the 
VMC method and constrained to have the correct asymptotic forms, produce S'-factor energy dependences for the 
processes ^H(a,7)^Li and '^He(Q;, 7)^Be which agree reasonably well with experiment. 

This work indicates no serious problems in the reaction rates presently used in astrophysical models. In fact, the 
most important implication of these results for astrophysics is probably that the previous understanding of these 
reactions remains essentially unchallenged. For example, using the present calculation to extrapolate the Robertson 
et al. activity measurement ||5^] of the S'-factor from 0.9 MeV to MeV, we obtain essentially the same result as 
with the energy dependence currently used in the standard solar model For big-bang nucleosynthesis, there is 
no low-energy extrapolation problem. The most useful result which a theoretical study could provide for cosmology 
would therefore be tighter constraints on cross-section normalizations than the current body of experimental data 
provides. The present calculations have not achieved that goal, but future first-principles calculations based on 
realistic nucleon-nucleon interactions might. 

The only serious problem with the results presented here is the low normalization of the ^'-factors, and its principal 
cause is probably easy to identify. Because most of the cross section arises at large (> 10 fm) cluster separations, the 
low normalizations most likely arise from the form of the seven-body bound states at large separations of the p-shell 
nucleons from the s-shell core. Part of the discrepancy may arise from the correlation fsp. The correlation between 
a and r clusters in the bound states is proportional to the twelfth power of fsp, so that the long-range correlation 
between clusters is very sensitive to the choice of fsp. However, it is hard to see how the fsp used in this study could 
affect the wave function beyond about 5 fm cluster separation. In the case of ^II(a, 7)^Li, a more likely explanation is 
that the Monte Carlo uncertainty has been underestimated, and a (prohibitively) long integration would "converge" 
on values in better agreement with experiment. Other aspects of the present calculation which are not done as exactly 
as one may wish, and therefore may bear on this problem, include the fact that asymptotic forms have been imposed 
on the wave functions which are inconsistent with their energies, as well as the use of cluster-cluster interactions which 
were phenomenologically constructed and could conceivably be inconsistent with the rest of the calculation in subtle 
ways. 

Despite these lingering difficulties, we have demonstrated the applicability of a new, almost ab initio approach for 
computing low-energy radiative captures in cases for which precise measurements exist, achieving the same accuracy 
as previous theoretical approaches. We also have presented additional evidence for uncertainty in the logarithmic 
derivative of the '^H(a,7)^Li S'-factor at zero energy. 

Regarding the present calculation as a first application of realistic potentials to radiative capture at mass seven, it 
has been successful, and it has taught important lessons about the role of the fsp correlations that were not apparent in 
previous VMC studies. This work clears the way for more refined models of these radiative captures based on realistic 
potentials. Specific improvements which will be possible in the near future include the use of improved three-body 
potentials now in development and the use of essentially exact wave functions derived by the Green's function 
Monte Carlo technique for the seven-body bound states. We note that the application of GFMC will require new 
wave functions, starting from VMC wave functions of the type used here to get the asymptotic forms correct. In the 
more distant future, it should be possible to perform the whole calculation self-consistently, constructing scattering 
wave functions from the potentials in a way similar to that used for the bound-state wave functions. VMC-based work 
would also profit by a more systematic effort to produce a Monte Carlo weighting scheme well-suited to computing 
the sorts of matrix elements encounted in direct capture calculations. Prospects for significant improvement on this 
initial investigation are very good. 
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